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Abstract 

We explore the effects of a heavy fermion doublet in a simplified version of the standard 
electroweak theory. We integrate out the doublet and compute the exact effective energy 
functional of spatially varying gauge and Higgs fields. We perform a variational search for 
a local minimum of the effective energy and do not find evidence for a soliton carrying 
the quantum numbers of the decoupled fermion doublet. The fermion vacuum polarization 
energy offsets the gain in binding energy previously argued to be sufficient to stabilize 
a fermionic soliton. The existence of such a soliton would have been a natural way to 
maintain anomaly cancellation at the level of the states. We also see that the sphaleron 
energy is significantly increased due to the quantum corrections of the heavy doublet. We 
find that when the doublet is slightly heavier than the quantum-corrected sphaleron, its 
decay is exponentially suppressed owing to a new barrier. This barrier exists only for an 
intermediate range of fermion masses, and a heavy enough doublet is indeed unstable. 
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I. INTRODUCTION 



In a chiral gauge theory like the Standard Model, gauge invariance prevents fermions from 
having explicit mass terms. Rather, they get their mass through their coupling to a scalar 
field via the well-known Higgs mechanism. At tree level the fermion obtains a perturbative 
mass, the product of the corresponding Yukawa coupling and the vacuum expectation value 
of the scalar Higgs field. The decoupling of such fermions presents interesting unsolved 
puzzles. Ordinary decoupling arguments which would show that a heavy fermion is 
irrelevant in the low-energy theory, break down. Increasing the mass, which causes the 
denominators in the fermion propagators to suppress quantum corrections, also increases 
the coupling, which gives a corresponding enhancement from the vertices. Furthermore, 
unlike an ordinary fermion with vector couplings, a chiral fermion cannot simply disappear 
from the theory as its mass is increased, because anomaly cancellation would be ruined. 
As shown in ref. gauge invariance is maintained at the level of the Lagrangian because 
integrating out the heavy fermion induces a Wess-Zumino term in the resulting effective 
Lagrangian. 

For the case of Witten's non-perturbative SU(2) anomaly 3], one can analyze the theory 
at the level of the action along the same lines M. However, one can also analyze the situation 

n 

from a different point of view. Ref. |4J shows that this anomaly can be understood in terms 
of the Hamiltonian of the theory: A theory with an odd number of left-handed fermion 
doublets has no gauge-invariant states. However, if the Yukawa coupling of a fermion is 
large enough, the perturbative fermion mass will be larger than the classical energy of the 
sphaleron 5J, so that such fermions are no longer stable states in the spectrum of the theory. 
Thus, to maintain gauge invariance in the low-energy theory, there must exist either new 
states in the theory carrying the quantum numbers of the fermion, or a mechanism to 
suppress the decay of the perturbative fermion states. 

Although these scenarios could rely on complicated non-perturbative physics, one simple 
resolution would be provided by the existence of a soliton carrying the quantum numbers of 
the decoupled fermion. If a localized configuration of gauge and Higgs fields binds a fermion 
level tightly, the binding energy could outweigh the cost in classical energy to set up the 
background field configuration. However, to consistently include the effects of the fermion 
level, such calculations must also include the Casimir energy, the renormalized shift in the 
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zero-point energies of all the other fermion modes, since both appear at the same order in 
h. For a static field configuration, the Casimir energy represents the full one-loop quantum 
vacuum polarization energy, equivalent to summing to all orders in the derivative expansion. 

To compare the quantum energy of the configuration with the sphaleron, we must also 
include the corresponding correction to the sphaleron's energy as well. Furthermore, we must 
check whether the quantum corrections induce an energy barrier between the perturbative 
fermion and the sphaleron, which would mean that the perturbative fermion would be quasi- 
stable, only able to decay by tunneling. 

We have carried out such calculations in a simplified version of the Standard Model, for 
background fields in the spherical ansatz, keeping fermion vacuum fluctuations but ignoring 
those of the gauge and Higgs fields. We find significant quantum corrections to the height 
of the sphaleron barrier. As we make the fermion level heavier than the quantum corrected 
sphaleron, we do see evidence for a barrier suppressing its decay. For even larger Yukawa 
couplings, however, the barrier disappears and the fermion's decay is unsuppressed. We 
do not see any evidence for a soliton for any value of the Yukawa coupling, and find that 
including the full Casimir energy destabilizes solitons found in previous work 

We comment on possible explanations of this result and its relation to the Witten anomaly 
in the conclusions. 

II. THE THEORY 

We consider the electroweak sector of the Standard Model with three simplifying mod- 
ifications: (1) the U(l) hypercharge gauge fields are absent, (2) the fermions within an 
isospin doublet are degenerate in mass, and (3) the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix is the identity, so there is no mixing between fermions of different generations. Since 
there is no U(l), the theory does not have any perturbative gauge anomalies, but Witten's 
non-perturbative anomaly requires it to have an even number of SU(2)l fermion doublets. 

The Higgs-Gauge sector Lagrangian density is 

C H = ~\tr (FTp,) + ^tr - ± [tr ($t$) _ 2v 2 ] \ (1) 

where 

Ffj,, = d^A v - d u A^ - ig [A^ A v ] , 
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= (d^-igA^)®, 

A* = ^y> ( 2 ) 

and where g and A are gauge and Higgs self-interaction coupling constants respectively. 
Furthermore, v denotes the tree-level vacuum expectation value of the Higgs field. The 
2x2 matrix field $ is related to the Higgs doublet <f> by 

(3) 

\ -01 02 J 

<3> can be re-expressed in terms of four real functions as 

$(ar) = v (s(x) + ip a (x)T a ) . (4) 

The fermionic sector Lagrangian density for each isospin doublet is 

C F = VLirfD^L + VRi-fd^R - f (V L $V R + h.c.) , (5) 

where the covariant derivative is the same as in eq. (J2J) and / is the Yukawa coupling constant 
(which may be different for each doublet). There is no coupling between fermions belonging 
to different isospin doublets owing to our assumption of diagonal CKM matrix. We introduce 
the potential 

V(A, $) = -grfA^x)!-^ + / (h(x) + ivp\x)T a l5 ) , (6) 

where 

h(x) = v{s(x) - 1) , (7) 
to write this Lagrangian density as the sum of a free part and an interaction part: 

C F = V [i^d^ -/«)*- Wtf . (8) 

We are interested in decoupling a single fermion doublet from the theory and for the re- 
mainder of the paper we consider only this single doublet. The full theory is defined by 

C = C H + C F . (9) 

In the unitary gauge (p a (x) = 0) at tree level, we have a single Higgs particle, h(x), 
with perturbative mass = 2v\f\. The superscript '(0)' denotes that the mass is at tree 
level. The three gauge fields are degenerate with perturbative mass = gv/ \/2. The two 
degenerate fermions have perturbative mass = fv. 
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III. THE HIGGS-GAUGE SECTOR EFFECTIVE ENERGY 



We integrate out the fermion doublet from the theory to obtain the effective action for 
the Higgs-Gauge sector: 

J[d^][d^}e i I dixC ^=o ■ 1 ; 

The normalization has been chosen so that the effective action is equal to the classical action 
for vanishing interaction potential, V, defined in eq. (JBJ). If iS^ denotes the Feynman 
diagram with one fermion loop and n external insertions of (—iV(A, $)), then 

/oo 
d 4 xC H + Y,SFD- (11) 
71=1 

The Feynman diagrams can be computed in a prescribed regularization scheme. The diver- 
gences that emerge as the regulator is removed are canceled by counterterms. We introduce 
the renormalized and counterterm Lagrangians, C^ en ^ and Cffl respectively, by expressing 
the bare parameters in in terms of renormalized parameters and counterterm coefficients 
to obtain Ch = C^ cn ^ + cffi. The renormalized Lagrangian is the original Higgs-Gauge sec- 
tor Lagrangian, eq. (0), with renormalized parameters substituted and 

C% t] = citr (F^F^) + c 2 tr 

+ c 3 [tr ($t$) _ 2v 2 ] + c 4 [tr _ 2v 2 ] 2 . (12) 

The coefficients c, depend on the regulator. For notational simplicity we do not introduce a 
different notation for the renormalized parameters. 

If we consider static Higgs-Gauge fields and restrict time to the interval T, then the 
effective energy functional is 

E cS [A, $] = - lim 1 S cS [A, $] = E cl [A, $] + E^A, $] , (13) 

T^OO _/ 

where E c \ refers to the classical energy of the Higgs-Gauge sector. We use the freedom to 
make time-dependent gauge transformations to set Aq = 0, so that 
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(14) 



E cl [A, $] = j d 3 x j~tr (FijFij) + itr ([A$] t D& ) + - [tr ($t$) _ 2w 2 ] 
The fermionic vacuum polarization energy is 

oo 

E vac [A, $] = ^ $] + ^ [A, $] , (15) 



71=1 



where each regulated Feynman diagram contribution is 

E&[4*] = - lim ±S<&, (16) 

1 —>oo 1 

and the counterterm contribution is 

B<* [A*] = y << 3 l|-Citr(FyFy) + c 2 tr 

-c 3 [tr ($ f $) - 2t> 2 ] - c 4 [tr ($ f $) - 2t> 2 ] 2 1 . (17) 

The entire one-loop effective energy receives contributions also from gauge and Higgs loops. 
We ignore these contributions because we believe the fermion loops are fundamentally re- 
sponsible for the phenomena associated with fermion decoupling. If we imagine that the 
fermions have N c internal degrees of freedom (e.g. color), then this approximation becomes 
exact for large Nc- Nevertheless we set N c = 1 in the analysis that follows. 



A. The Counterterms 

The counterterms render E vac finite by canceling the divergences in E^, for n — 1 
through n — 4, that arise when the regulator is removed. To unambiguously determine the 
counterterm coefficients, q, we impose conventional renormalization conditions: 

a. We choose the vacuum expectation value (vev) of h(x) to be 0, which ensures that 
the vev (<£>) = vl stays fixed at its classical value and the perturbative fermion mass 
does not get renormalized, i.e. rrif — m^. This is equivalent to a "no-tadpole" 
renormalization condition and determines c 3 . 

b. We fix the pole of the Higgs propagator to be at the tree level mass, m h = mf \ with 
residue 1. These conditions yield c 2 and c 4 . 

c. We have various choices to fix the remaining undetermined counterterm coefficient c\. 
We choose to set the residue of the pole of the gauge field propagator to 1 in unitary 
gauge. Then the position of that pole, i.e. the mass of the gauge field, m w , is a 
prediction. 

The resulting counterterm coefficients, q, are listed in Appendix A. As explained under item 
c, the mass of gauge fields is constrained by the other model parameters when fermion loops 
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are included. With our choice of renormalization conditions, it is the solution to the implicit 
equation 



K 0) ) 
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dxx (1 — x) 



A(x,ml 1 ) 



i 1 , , m A(x,m 2 h ) 
+6 / dxx(l-x)hy- v " 
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with A(x,g 2 ) = m 2 — x(l — x)g 2 . Recall that mL = gv/y/2 is the tree level perturbative 
mass of the gauge fields. 



B. The Vacuum Polarization Energy 

We briefly summarize methods introduced in refs. 13, Q (see ref. for a review and 
a list of additional references) that enable an exact calculation of the fermionic vacuum 
polarization energy. Eq. (|T5j) gives this energy as an infinite sum of Feynman diagrams. 
For non-perturbative field configurations all orders would have to be summed, which is 
intractable. We therefore make use of the fact that it is also given by a sum over the shift 
in the zero-point energies of the fermion modes due to the background fields. We write this 
formal quantity as a sum over bound state energies, ej, (times their degeneracies, Dj) and 
a momentum integral of the energy of the continuum states weighted by the change in the 
density of continuum states, Ap(k), that is induced by the background fields, 



E, 



vac 
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-J^D^l--^ dk^k* + m}Ap(k)+E ct . (19) 



The momentum integral and E ct are both divergent, but their sum is finite because the 
theory is renormalizable. We render the integral finite by subtracting the first N terms 
in the Born series expansion of the density of states and adding back in exactly the same 
quantity in the form of Feynman diagrams: 



N 

+ Y^E^ + E ct . (20) 

i=l 

When we combine E ct with Y2iLi -^fem we cancel the divergences as well as implement the 
renormalization prescription. As a result, the above expression is manifestly finite. The 



minimal number of required Born subtractions, N, is easily determined by an analysis of 
the superficial degree of divergence of the one-loop Feynman diagrams. For our theory, 
N = 4. 

I — I 

We will work with background fields in the spherical ansatz |lOj ]. Then we can express 
Ap(&) (and its Born series) in terms of the momentum derivative of the phase shifts jll|. 
induced by the background fields, of the Dirac wave-functions, 

<re{+,-} G 

Here we have expanded the 5-matrix in grand spin channels labeled by G. The grand spin 
is the vector sum of total angular momentum and isospin. It is conserved by the potential, 
eq. ©• The asymptotic scattering states are labeled by parity (— 1) G and total spin G± 1/2 
and we obtain a 4 x 4 S- matrix in general (except in the G = channel, where it is 2 x 2). 
We let Sa,a{k) denote the sum of the eigenphase shifts at momentum k in channel G and 
a — ± specify the sign of the energy eigenvalue, uj = ±^/A; 2 + raj. Note that the single 
particle spectrum is not symmetric because the Dirac Hamiltonian is not charge conjugation 
invariant. The degeneracy is given by Dq = 2G + 1. In Appendix B we show in detail how 
to use the Dirac equation to calculate the bound state energies and scattering phase shifts 
needed in the computation of E v3bC . 

To simplify the calculation we only subtract the first two Born approximants and elim- 
inate the remaining log- divergence in the momentum integral by using a limiting function 
approach developed in ref. [12j. The final expression for the vacuum polarization energy is 

£vac = ~ £(2G, + l)|e,| -[^ Jk* + m*±6(k) + + E^ , (22) 
where Gj is the grand spin associated with the bound state j and 

oo 

E E( 2G+1 )(^(*)-C(*)-C(*))+w*)- (23) 

ae{+-}G=o 

Here S^ a (k) denotes the z th -term in the Born series of 8a,a{k)- After subtracting <$W and 
5(2) f rom 5 ) the momentum integral does not contain any contributions that are linear 
or quadratic in V. The limiting function for the sum over all eigenphase shifts, 5u m (k), 
eliminates the logarithmically divergent pieces that are third and fourth order in V from the 
momentum integral. Its analytic expression can be extracted from the divergent pieces of 
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the corresponding Feynman diagrams and is given in Appendix B, eq. (|B18|) . Furthermore, 
E^ 1,2 ) denotes the contribution up to second order in V from the renormalized Feynman 
diagrams. Its explicit expression is displayed in Appendix B. Finally, E^ 3 ' 4 ^ contains the 
third and fourth order counterterm contribution combined with the divergences in the third 
and fourth order Feynman diagrams that have been subtracted from the momentum integral 
via S\i m . Again, its explicit form can be found in Appendix B. 

Thus we compute an exact, finite, renormalized, gauge-invariant effective energy func- 
tional, E e g[A, <&], defined in eq. (fT^j). for the Higgs-Gauge sector, with the fermion fields 
integrated out. 

IV. THE ENERGY OF A FERMIONIC CONFIGURATION 

We are interested in exploring the possibility of the emergence of a stable, fermionic 
soliton in the Higgs-Gauge sector of the theory as we increase the Yukawa coupling (thereby 
making the perturbative fermion heavier). In the previous section we outlined the procedure 
that allows us to calculate the effective energy when the fermions are integrated out. Now 
we analyze the minimum additional energy required to associate unit fermion number with 
a particular Higgs-Gauge configuration C, where C = {A, $}. First in section Ft V Al we 
determine the integer fermion number F[C\ of a configuration. (This is subtle because we 
have to contend with the anomalous violation of fermion number). Then we occupy or 
empty levels of the single-particle Dirac Hamiltonian to obtain the lowest energy state with 
net fermion number 1. If F[C] — 0, then the lightest positive bound state needs to be filled 
and the occupation energy Eoll = ei owe st, where the superscript '(1)' denotes that levels have 
been occupied/emptied to obtain fermion number 1. If F[C] = 1 then Eoll = because C is 
already fermionic, and so on. Thus, the minimum total energy of a single fermion associated 
with a configuration is 

[C] = E cl [C] + £ vac [C] + E<» [C] . (24) 

In previous works, such as jf|, the vacuum polarization contribution was omitted from the 
above equation and stable non-topological solitons were found. We consider such calculations 
inconsistent, because Eoll and E vax . are both order h. We will see explicitly in section IV~Gl 
that E v3iC makes a significant positive contribution when the Yukawa coupling is large enough 
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that the perturbative fermion mass is comparable to the classical energy. 

Since we refer to all these different energies frequently in the rest of the paper, we sum- 
marize our notation in Table 





Classical Higgs and gauge energy, eq. (fTl)) 


E 


Fermion vacuum polarization energy, including counterterms, eqs. (|15l I19|) 


E e S 


One-fermion-loop effective energy, E c \ + E va/C 


P ( m ) 
-G/occ 


Smallest occupation energy required to obtain fermion number to, eq. 1)24(1 


E C S 


Smallest effective energy in the fermion number to sector, Egg + Eo™c 



TABLE I: Definitions of some of the energies which appear in our analysis. 



A. The Fermion Number of a Configuration 

First we review properties of the Higgs- Gauge configuration space and the classical energy 
functional defined on it. From the expression for E c \ in eq. (|14|) . it follows that configurations 

A, = -U^djU^ , $ = vU {n) , (25) 
9 

have E c \ = and we refer to them as zero- classical-energy configurations. Here [/(") is any 
map from S 3 to SU(2) with winding number n. Zero-classical-energy configurations with 
the same winding number are equivalent under small (winding number 0) gauge transfor- 
mations. We use to denote the homotopic class of zero-classical-energy configurations 
with winding number n. Topologically inequivalent zero-classical-energy configurations are 
related by large (nonzero winding number) gauge transformations. Along any continuous 
interpolation between two configurations, one in and the other in (with n ^ to), 
there exists a configuration, C, with maximum classical energy. Since no can be contin- 
uously deformed into any U^ m \ E C \[C] > 0. The configuration corresponding to the minimax 
of these energies, when all interpolations are considered, is the classical sphaleron. When the 
fermion vacuum polarization energy is added to the classical energy to obtain the effective 
energy (E eS = E c \ + E VSLC ), not only does the magnitude of the minimax energy change, but 
its location in configuration space shifts as well. We therefore define the quantum- corrected 
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sphaleron to be the configuration that has the lowest of the maximum effective energies along 
all interpolations between topologically inequivalent zero-classical-energy configurations. 

We associate any configuration C with a unique class of zero-classical-energy configura- 
tions by continuously deforming C in the direction of the negative gradient of the classical 
energy until we get a configuration in for some n = n(C). We call C^ c '' the connected 
C- class of C and we say that C is in the winding number n basin. Note that the classical 
sphaleron and all configurations that descend to it are on the boundary between different 
basins. 

For any two configurations G\ and C2, we define the spectral flow S[Ci,Cq\ to be the 
number of eigenvalues (levels) of the single particle Dirac Hamiltonian that cross zero from 
above minus the number that cross zero from below along any interpolation from C\ to 0%. 
The fermion number of a configuration C is defined as 

F[C]=S[C 1 ,C] (26) 

with C\ E C^ c ''. Since the Dirac spectrum is gauge- invariant, F[C] does not depend 
on which particular C\ is chosen from the connected C-class. Moreover, F[C] is gauge- 
invariant even under large gauge transformations. Also, F[C] does not depend on the chosen 
interpolation because the spectral flow is the same for all interpolations. This definition of 
the fermion number can be readily understood for a C that has classical energy less than 
the energy of the classical sphaleron. A continuous interpolation from any configuration in 
C^ c » to C preserves net fermion number because anomalous fermion number violations 
require the Higgs-Gauge fields to cross the sphaleron barrier. Thus, defining zero-classical- 
energy configurations to have zero fermion number leads to eq. (]26[). Configurations that 
have classical energies larger than the classical sphaleron are not separated by an energy 
barrier from topologically inequivalent basins, so it is not clear what their fermion number 
should be, although our definition does assign a unique F to them. 

Having determined the fermion number of a configuration, we can use eq. (J2%j) to find 
, the minimum effective energy of a configuration in the fermion number 1 sector. 



11 



B. Stability of the Soliton 



We would like to know if there exists a configuration at which the one-fermion effective 
energy functional, , has a local minimum. This configuration would be a fermionic soli- 
ton. We carry out a variational search, looking for a configuration C such that E^[C] < nif 
and £gjj [C] < E q , s _, where E q _ s , is the effective energy of the quantum-corrected sphaleron. 
The first condition ensures that C cannot simply decay into a configuration with zero clas- 
sical energy plus a perturbative fermion. The second condition prevents C from rolling 
over the quantum-corrected sphaleron, giving up its fermion number and then rolling down 
the E^ surface to a zero-classical-energy configuration. Finding a configuration with these 
properties would guarantee the existence of a nontrivial local minimum of E^ . 



V. THE SEARCH FOR THE SOLITON 



In this section we describe our search for the soliton. We first review the spherical ansatz 
for the gauge and Higgs fields. We then outline the restrictions imposed on the variational 
ansatze used to search for a soliton. Finally, we report on our search within two physically 
motivated sets of ansatze: the "twisted Higgs" and "paths over the sphaleron" . Note that 
throughout this section the perturbative fermion mass is set to 1 so that energies and lengths 
are measured in units of mj and 1/rrif, respectively. 



A. The Spherical Ansatz 



We only consider static gauge and Higgs fields in the spherical ansatz. This enables us 
to expand the fermion S matrix in terms of partial waves labeled by the grand spin G, 
as described in Appendix B. Our method for calculating the fermion vacuum polarization 
energy requires such an expansion. Under these restrictions (and in the A = gauge, 
which for smooth fields is obtained by a non-singular gauge transformation), the fields can 
be expressed in terms of five real functions of r: 

1 



$(£) 



-A l (x) 



/ s „ . a(r) . „ „ . 7(r) 



v [s(r) + ip(r)rjXj] , 



(27) 



where x is the unit three-vector in the radial direction. 
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The ansatz transforms under a U(l) subgroup of the full SU(2) gauge symmetry, corre- 
sponding to elements of the form 



9{x) 



(28) 



with ai transforming as a 1 dimensional vector field, s + ip as a scalar with charge 1/2, and 
a + 2(7 — 1) as a scalar with charge 1. It is thus convenient to introduce the moduli p, S 
and phases 9, 77 for the charged scalars: 

- ipe ie = a + i{j - 1) and Ee iv = s + ip. (29) 

From now on we will specify a configuration using the five functions ai(r), p(r), 9(r), £(r) 
and ?](r). Regularity of and at x = requires that 

P(0) = 1, 

p(o) = 0, 

9(0) = -2n e iT, 
ai(0) = 0'(O), 

either E(0) = or 77(0) = -n^TT. (30) 

Here ng, n v are integers and primes denote derivatives with respect to the radial coordinate. 

For the gauge transformation g(x) in eq. (|2*8*|) to be non-singular, we require f(0) = —2nir, 
where n is an integer and we denote this boundary condition as a superscript: /(r) = f^ n \r). 
If f( n '(r) is restricted to be as r — > 00 (which is equivalent to g(r — > 00) = 1) then n is 
the winding of the map g(of) : S* 3 — ► SU(2). So the topology of the zero-classical-energy 
configurations persists in the spherical ansatz. The classical energy of eq. (j!4|) takes the 
form 



E, 



cl 



47T 



19 



p 2 (9'-a 1 ) 2 + 



(P 2 " 1) 



2r 2 



+ 



./ 2 



r 
1 



(p-l) 2 + 4p 2 sin 2 ^ 2V 



2 V 2 
and winding number n zero-classical-energy configurations of eq. (|25jl now become 



(31) 



p(r) 



1 



= / (n) (r), 7?(r) 



/( n) (r) 



«i(r) = / (n)/ (r). 



(32) 
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We want the Higgs-Gauge fields to have finite classical energy. So we require a field 
configuration of the form eq. (|32|) as r — > oo, and the restriction that /^(oo) = uniquely 
specifies the boundary conditions on the fields at infinity. At x — 0, the boundary conditions 
on p (specified in eq. (|3UJl ) make the energy density finite and we do not require any additional 
constraints. 

The anomalous violation of fermion number is given by the anomaly equation 

d,(^^)=d,K», (33) 

where is the Chern-Simons current. It is useful to consider the charge associated with 
it: 

^cs = - J^'fc J d3xtr (y A i d i A k ~ ^igAiAjA^j 

= — / dr [a 1 + p 2 (9'-a 1 )] . (34) 
Z7[ Jo 

This is a non- integer in general, and is equal to the integer winding number of for 



the configurations of eq. (|32|). and a half- integer for the sphaleron 13] . Under a winding 
n gauge transformation, iVcs — ► Ncs + n. For background fields that interpolate between 
topologically distinct zero-classical-energy configurations, the net fermion number produced 
is given by the change in Nq§. 



B. Restrictions on the Variational Ansatze 



Our methods allow us to consider any static, spherically symmetric configuration, C, in 
the Higgs-Gauge sector, specified by the five real functions aj.(r), p(r), 0(r), E(r) and rj(r). 
In principle, we could numerically minimize the fermionic energy, E^[C], in terms of the 
five functions and determine if a soliton exists. In practice however, that is not feasible. 
So instead we limit ourselves to the variation of a few parameters in ansatze motivated by 
physical considerations. 

We restrict our variational ansatze to those that obey the above described boundary 
conditions at the origin and at infinity. In addition, we restrict the Higgs fields to lie within 
the chiral circle, E(r) < 1, because otherwise the effective potential is unbounded from 
below. (The leading terms in the derivative expansion of eq (10) can be found in Q|). 
Finally, the effective theory (obtained by integrating out the fermions) has a Landau pole 
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in the ultraviolet, reflecting new dynamics at some cutoff energy scale or equivalently at 
a minimum distance scale. Configurations that are large compared to this distance scale 
are relatively insensitive to the new dynamics at the cutoff, but smaller configurations are 
sensitive. For small widths and large couplings, the Landau pole becomes significant and 
leads to unphysical negative effective energies in eg. (|T5|l is1 |. We have to be wary of this 



in estimating the reliability of our results. See 
effective potential and the Landau pole. 



121 for more detailed discussions on the 



C. Twisted Higgs 

We first consider twisted Higgs configurations, with n r] = 1 so that rj(r) goes from — it 
at r = to as r — > oo. The other functions are trivial: a\(r) = 0,p(r) = l,8(r) = 
and S(r) = 1. If we smoothly interpolate from a zero-classical-energy configuration (in the 
connected C-class) to such a configuration, we observe that one fermion bound state, that 
originates in the positive continuum, has its energy decrease sharply. The wider the final 
twisted Higgs configuration, the closer this level ends up to the negative continuum. At a 
width of order l/rrif, it has energy zero, eliminating any occupation energy contribution, 
-Eocc, to the fermionic energy, E^g , associated with it. The existence of this level makes such 
twisted Higgs configurations attractive candidates for the variational search. 

We consider one such twisted Higgs configuration with a width characterized by a varia- 
tional parameter w, 

v = -ne~ r/w (35) 

and add various perturbations to it. For instance, one among many of our variational ansatze 
(in the 8 = gauge) is a four parameter ansatz with parameters Po, . . ■ ,ps. 

-r/w i r / w -r/w 



1 + (r/w) 

£ = 1+ ^ T-7 / '^ r/w > 

1 + (r/w) 

r / W -r/w 

1 + (r/w) z 
1 + (r/w) 6 

where —1 < pi < (to keep the Higgs field within the chiral circle and its magnitude 
positive) and P3 > —5.23 (to keep p positive). For a prescribed set of theory parameters 
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(m w ,mh and /) we determine the gauge coupling g = ^/2mw /v from the renormalization 
constraint eq. (|TH)l. We then vary the ansatz parameters (w, p , . . . , p 3 ) to lower the fermionic 
energy . We find that the gain in binding energy is insufficient to compensate for the 
increase in the effective energy E e ^. Through all our variations, is strictly greater than 
rrif and we find no evidence for a soliton. The same result was obtained in |l2[ without 
gauge fields, and the extra gauge degrees of freedom do not seem to help in the twisted 
Higgs ansatz. 

We find that the fermion vacuum polarization contribution, E vac , to the fermionic energy, 
E^ destabilizes would-be solitons. Consider a linear interpolation from the trivial zero- 
classical-energy configuration to the twisted Higgs configuration in eq. (J35|) with gauge fields 
set to zero. We introduce the interpolating parameter £ which goes from to 1: 

Ze iv = 1 - £ + £ exp (-tire- r/w ) . (37) 

We choose the Yukawa coupling to be / = 10 and the Higgs mass to be v/y2. Since the 
gauge fields are trivial, the classical energy as well as the Dirac spectrum are independent 
of the gauge coupling g and the gauge bosons mass m w . For each value of £, we compute 
Etff and E^ — E vac for different values of the width parameter w. In Fig. ^we plot E^ 
and E^ — E vac as functions of £, choosing the width at every point to minimize the energy 
(we do not allow w to be less than 1 so that we remain relatively insensitive to the Landau 
pole). For all points on the plot there is no spectral flow and so E^ = E c \ + e lowcst + E yaic in 
accordance with eq. (j2U), where ei OW est is the smallest positive bound-state energy in the Dirac 
spectrum. If E vac is omitted, for < p < 0.6 we have configurations that have fermionic 
energies lower than the mass of the perturbative fermion (m/ = 1 in our units). These 
configurations indicate the existence of a local minimum on the E^g — E vac surface which 
would be a soliton. The E vac contribution, however, raises the energies of the configurations 
to above m/ as shown in the figure, and the would-be solitons are destabilized. 

D. Paths over the Sphaleron 

The gauge fields introduce another mechanism for strongly binding a fermion state: there 
is a zero mode in the background of the sphaleron. Along an interpolation of the background 
fields from a configuration in to a configuration in C^ 1 -*, a fermion level leaves the 
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FIG. 1: Minimum fermionic effective energies (in units of mf), with as well as without E vac 
contributions, along the interpolation in eq. (|37j) . 

positive continuum, crosses zero from above and finally enters the negative continuum. The 
lowering of the occupation energy, Eoll, as we approach zero from above is offset by the rising 
effective energy (E e R = E c \ + E vac ), so we must investigate whether the former can dominate 
the latter. We also use such interpolations to study the effects of a large Yukawa coupling on 
the sphaleron. We approximate the quantum-corrected sphaleron by minimizing the effective 
energy barrier between topologically inequivalent zero-classical-energy configurations, with 
respect to the variational parameters of our interpolations. We also investigate the possible 
emergence of new barriers in the one-fermion energy surface when the perturbative fermion 
becomes heavier than the quantum-corrected sphaleron. These last two phenomena affect 
the stability of the heavy fermion and in some models may be significant for baryogenesis. 

We make the following choices for the theory parameters: we fix the Yukawa coupling 
at / = 10, which is large enough that fermion effects are significant, but small enough to 
prevent our configurations from being affected by the Landau pole. Indeed we encounter 
no negative energy instabilities in our computations of _E vac for this coupling. We keep the 
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Higgs mass fixed at v/y2, which corresponds to rrih = -^j ~ 0.07m/. We choose g to keep 
the quantum-corrected sphaleron energy comparable to m/, since it is plausible that the 
mass of the sphaleron sets the scale for decoupling. If the sphaleron is too heavy compared 
to the fermion, the binding energy gained by the level crossing would be washed out by 
the effective energy of the sphaleron. For g = 6.5 our best approximation to the quantum- 
corrected sphaleron is approximately degenerate with the fermion. When the gauge coupling 
is given, the mass m w of the gauge boson is determined from the renormalization constraint 
eq. ()18|): m w ph 0.63m f for g = 6.5 and m w « 0.98m f for g = 10. These theory parameters 
are of course large deviations from the Standard Model parameters. We exaggerate them 
to see the effects of the heavy perturbative fermion. Another concern is that for large g, we 
should consider quantum fluctuations of the Higgs-Gauge fields. However, we believe that 
anomaly cancellation drives the creation of a soliton, which would suggest that the fermion 
vacuum polarization contains the essential physics, and our methods allow us to exactly 
compute this contribution to the energy for any Yukawa coupling. 

First we consider a linear interpolation between a winding-0 and a winding- 1 zero- 
classical-energy configuration. The interpolation parameter £ goes from to 1: 



In the spherical ansatz, U^(x) = g(x) is specified by a single function, as in eq. (|2*8*|) . that 
we choose to be 



where w characterizes the width of the configuration. We interpolate £ from (the trivial 
configuration with Acs = 0) to 1/2 (a configuration with Acs = 1/2 in the presence of which 
the fermion has a zero mode) and vary w along the interpolation to minimize E^ . This 
gives an upper bound on the minimum as a function of Acs- For Acs = 1/2, this is an 
upper bound on the quantum-corrected sphaleron energy as well, because E^g = E eS in the 
presence of a fermion zero mode (the occupation energy is then 0) 2 . An exploration of Acs 

2 Our numerical methods do not allow us to consider Ncs exactly equal to 1/2, because the Higgs magnitude 
vanishes at r — and the second order Dirac equations develop a singularity (as discussed in Appendix 
B), but we can compute E^g very close to Ncs = 1/2- 



$ = v(l-£)l + £vU m , 
A j = tiuWdjUW* . 



9 



(38) 




(39) 
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between and 1/2 is sufficient to map to all values of Nqs, since configurations with Ncs 
between 1/2 and 1 are obtained by charge conjugation, and configurations with iVcs < or 
iV cs > 1 are large-gauge-equivalent to configurations with N cs between and 1. 

We also consider an instanton-like configuration where the Euclidean time £ = 24 is the 
interpolation parameter (which varies from —00 to 00) between two topologically inequiva- 
lent zero-classical-energy configurations: 

A, = h{r,0~U in 4x,t)d»Ul st {x,0, 
$ = vVh(r,£)U inBt (x,£), (40) 

where 

U^(x,0 = ^^ 2 (41) 
is the canonical winding-1 map from S 3 (space-time infinity) to SU(2). Furthermore h is a 



function of the Euclidean space-time radius (y/r 2 + 
goes from to 00. 't Hooft's electroweak instanton 



and goes from to 1 as this radius 
is constructed clS db self-dual gauge 



field configuration in the topological charge one sector, and a Higgs field configuration that 
minimize the covariant kinetic term in the Lagrangian density. This gives 

r 2 + £ 2 

h <r>Z)= 2IT2I 2 ( 42 ) 

for any width w (the classical theory with no Higgs field is scale- invariant). We modify this 
radial function to exponentially approach its asymptotic value of 1, so that the potential 
in our Dirac equation falls off fast enough to have a well-defined scattering problem (as 
described in Appendix B). We choose 

h(r,£) = 1 - e -( r2+ « 2 ^ 2 . (43) 

This choice does not minimize any part of the classical Euclidean action (in the topological 
charge one sector). Since we are interested in minimizing , which has fermion vacuum 
energy and occupation energy contributions in addition to the Higgs-Gauge sector classical 
energy, we do not need our configurations to minimize the classical energy. In fact, as we 
describe later, the configurations that minimize are rather different from those that 
minimize E c \. 

In order to compute the Dirac spectrum for this background using the methods described 
in Appendix B, we gauge transform to A = and lim^oo r\ = lim^oo 9 = using the 
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FIG. 2: Minimum effective energies (in units of m/) along the linear path in eqs. (j381 1.39|) . in both 
the zero-fermion and one-fermion sectors (with as well as without the E vac contributions). 



transformation function 



f(r,0 



ft 2r 



(44) 



in eq. (|28|). Finally, as in the case of the linear interpolation, we consider £ from — oo (the 
trivial configuration with Nqs = 0) to (a configuration with iVcs = 1/2 in the presence of 
which the fermion has a zero mode) and for each £ we choose the w that minimizes E^ . 

In Fig. 121 we plot the minimum effective energies in both the zero-fermion sector (E^) 
and the one-fermion sector (E^), minimized within our variational ansatz for the linear in- 
terpolation, as functions of the Chern-Simons number, iV cs (see eqs. (I38I39J) ). As mentioned 
before, we fix the theory parameters at a Yukawa coupling of / = 10, a gauge coupling of 
g = 6.5 and a Higgs mass of ~ 0.07mj. These parameters determine the mass of the 
gauge bosons to be m w w 0.63m f from the renormalization constraint, eq. (|18|) . To isolate 
and highlight the contribution of the fermion vacuum polarization energy, E vac , we also plot 
the energies minimized with E vac subtracted. 

For 



First consider the zero-fermion sector with the two curves E^ and E^' — E y 
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all points on the plot there is no spectral flow and so — E v3iC = E c \. At Ncs = 0, 
both E^g and E c \ are minimized at the trivial zero-classical-energy configuration, eq. (|32|) 
with f^(r) = 0. At iVcs = 1/2, E^ is minimized at the quantum-corrected sphaleron 
while E c \ is minimized at the classical sphaleron. Within our variational ansatz, we find the 
parameters that minimize E^ at Ncs = 1/2 are different from those that minimize E c \. So 
our approximation to the quantum-corrected sphaleron is distinct from our approximation 
to the classical sphaleron. Moreover, the fermion vacuum polarization energy correction 
to the sphaleron turns out to be rather large. Our classical sphaleron has an energy of 
0.45m/ (which agrees well with the numerical estimate of E = 1.52 47TV ^ in [13]), while our 
quantum-corrected sphaleron has an energy of 1.02m/. 

Next consider the E^g and E^g — E vac plots in the one- fermion sector in Fig. El Again, for 
all points on the plot there is no spectral flow and so E^ = E c \ + ei owost + E vac in accordance 
with eq. (J2U), where ei owest is the smallest positive bound-state energy in the Dirac spectrum. 
Since the classical sphaleron has an energy much smaller than the perturbative fermion mass, 
one would expect that the perturbative fermion would have an unsuppressed decay mode 
over the sphaleron, as first pointed out by Rubakov in 3). The — E vac curve indeed 
displays this decay path. The fermion vacuum polarization energy modifies things in two 
crucial ways. First, the fermion quantum corrections to the sphaleron raise its energy to 
be degenerate with the fermion, as mentioned before. So the threshold mass is significantly 
increased. Second, in the plot of E^ we observe that there is an energy barrier between 
the fundamental fermion and the quantum-corrected sphaleron. This indicates that even 
when the fermion becomes heavier than the sphaleron, there might exist a range of masses 
for which the decay continues to be exponentially suppressed (since it can proceed only via 
tunneling) . 

In Fig. El we restrict our attention to the one-fermion sector and consider the linear 
interpolation and the instanton-like interpolation. In addition to g = 6.5 we consider g = 10, 
which corresponds to m w m 0.98m f. For each of these two gauge couplings, we plot the 
effective energy minimized within our ansatze as a function of Nqs for the two interpolations. 

The two seemingly different interpolating configurations produce very similar minimum 
E^q curves. Furthermore, when we enlarge the variational ansatze in our interpolations, 
we are unable to reduce the energies by any significant amount. Thus we speculate that 
the plotted curves may be close to tight upper bounds on the true minimum E^ curve. 
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FIG. 3: Minimum (in units of mf) for paths from a zero-classical-energy configuration to the 
sphaleron. Curves denoted 'Lin.' refer to the linear path in eqs. (|3*51l39|) while those labeled 'Inst.' 
are associated with the instanton path, eqs. (|4UI41I43|) , 

This is the justification for considering only the linear interpolation in Fig. |2l and taking 
the evidence for the emergence of a new barrier and the significant energy change of the 
sphaleron seriously. 

Note that as the gauge coupling increases, lowering the energy of the quantum-corrected 
sphaleron, the barrier between the fundamental fermion and the sphaleron does not persist 
indefinitely. We observe that for g = 10, when rrif is approximately 1.3 times our quantum- 
corrected sphaleron, there is no barrier and the decay mode is finally unsuppressed. 

Just as in the case of the twisted-Higgs variational ansatz, in our ansatze of paths from a 
zero-classical-energy configuration to the sphaleron, we have not found a configuration with 
the associated fermion energy lower than both the perturbative fermion and the quantum- 
corrected sphaleron. Thus, we find no evidence for the existence of fermionic solitons in the 
low energy spectrum of the Standard Model within our ansatze. 



22 



VI. CONCLUSIONS AND DISCUSSION 



We have explored quantum effects of a heavy fermion on a chiral gauge theory. The 
quantum-corrected sphaleron is heavier than the classical sphaleron by an energy of the order 
of the perturbative fermion mass. This higher barrier suppresses fermion number violating 
processes. We also observe the creation of an energy barrier between the perturbative 
fermion and the sphaleron, so that a fermion with energy slightly above the sphaleron can 
still only decay through tunneling. This new barrier exists only for an intermediate range of 
perturbative fermion masses, and a heavy enough perturbative fermion is indeed unstable. 
We do not, however, find evidence for the existence of a soliton in the spectrum of the theory. 
The fermion vacuum polarization contribution seems to destabilize any would-be solitons. 

It is possible that for a large enough Yukawa coupling, the Witten anomaly is saturated 
by states that do not have a particle interpretation. But we believe the anomaly puzzle 
could still be resolved by soliton states that exist outside the spherical ansatz. Although the 
reduced U(l) theory of the spherical ansatz reproduces the anomalous violation of fermion 
number, it does not have the Witten anomaly. This suggests that the ansatz might be too 
restrictive to resolve decoupling issues associated with the Witten anomaly. Klinkhamer has 
conjectured the existence of a non-spherically-symmetric sphaleron jl^ derived using the 
same topological properties underlying the Witten anomaly, which would provide one can- 
didate background. Like the usual sphaleron, this configuration has a zero mode because it 
represents the middle of a path in which a fermion level crosses zero. However, an evaluation 
of the stability of such a state would require a calculation of the one-loop effective energy, 
which is more difficult in the absence of spherical symmetry. 

ACKNOWLEDGMENTS 

We would like to thank J. Goldstone, M. Quandt, and F. Wilczek for valuable discus- 
sions and also P. Sundberg for his help. E. F., R. L. J. and V. K. are supported in part by 
the U.S. Department of Energy (D.O.E.) under cooperative research agreement #DF-FC02- 
94ER40818. N. G. is supported by the U.S. Department of Energy (D.O.E.) under coop- 
erative research agreement #DE-FG03-91ER40662 and H.W. is supported by the Deutsche 
Forschungsgemeinschaft under contract We 1254/3-2. 



23 



APPENDIX A: RESULTS FROM FEYNMAN DIAGRAMS 



In this Appendix we list the results from the Feynman diagram calculations mentioned 
in section III. 

In dimensional regularization (d = 4 — e) the counterterm coefficients, defined in eq. (|12j). 
read 
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where /i is the momentum scale introduced to maintain the canonical dimensions of the 
parameters when regularizing in fractional dimensions. 

In eq. (I22|) E^- 1 ' 2 ' denotes the contribution to the vacuum polarization energy from first 
and second order renormalized Feynman diagrams. Its explicit expression reads 
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with the Fourier transform of a field <p(x) defined in the usual way as <p(q) = J d 3 xip(x)e tq ' x . 
The third and fourth order counterterm contribution combined with the divergences in the 
third and fourth order Feynman diagrams is 
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where = $ — t> parameterizes the deviation of the Higgs field from its vev. 



APPENDIX B: THE DIRAC EQUATION 

In this Appendix, we describe how we obtain the bound state energies of the Dirac equa- 
tion and the scattering phase shifts (and their Born series) in the presence of a background 
potential. These quantities are required to compute the vacuum polarization energy in 
eq. (ED- 

The fermion field obeys the time-independent Dirac equation 

H D ^ = , (Bl) 

where 

H D = -i-fydt + 7° K + V(A, $)] , (B2) 

and V is given in eq. (JSJ). In contrast to previous work, it is most convenient to use the 
chiral representation of the Dirac matrices, 

i<Jjdj + g<JjAj m^(s + ipTjXj, 
rrif(s — ipTjXj) —idjdj 
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The grand spin G is defined as the vector sum of isospin, spin and orbital angular momen- 
tum. It commutes with Hp as long as the fermion doublet is degenerate in mass and the 
background fields are in the spherical ansatz. We satisfy both conditions. For a given grand 
spin quantum number G (we suppress the grand spin projection label M throughout), the 
Dirac spinor ty G has eight components and may be written in terms of generalized spherical 
harmonic functions 3^,i(£) with j — G ± 5 and I = j ± | as 

= ( *fAi + 92y G+hG + g 3 y G - hG + ig4y G -\,G-i \ ? (B4) 
y ihy G+ ^ G+l + /2^ G+ | )G + hy G - i, G + ihy G -^ G -! ) 

where <?j and fi are radial functions and we have suppressed the grand spin labels on them. 
Note that in this chiral theory modes of different parity, e.g. g\ and g 2 mix. The spherical 
harmonics are two-component spinors in both spin and isospin space. The special case \l/o 
is defined only in terms of 34^ and and does not contain g 3 , g 4 , / 3 , / 4 . 

The matrix elements of operators like TjXj between the spherical harmonics may be 
found in the literature [3]. We use them to write the Dirac equation (|B2|) as a a set of eight 
coupled first-order linear differential equations in the radial functions, for fixed G. From 
these equations we obtain the bound state solutions (\uj\ < rrif) in each grand spin channel 
using shooting algorithms. From Levinson's theorem we determine the number, N G ound , of 
bound states to shoot for, using phase shifts, S G (cu), of the scattering state solutions of the 
Dirac equation: 

^bound = 1 (5c(m/) _ 5g(qo) + s G (- mf ) - 5 G (-oo)) . (B5) 

7T 

To construct these scattering state solutions we re-write the Dirac equation as a set of 
second-order differential equations in the radial functions. Formally they read, 

[h 12 h 21 - h l2 {h 22 - u)h${h u - u;)] ttg = , (B6) 

where 

= my G+ i, G +i + 92y G +\,G + 9sy G -i, G + my G -±, G -i 

denotes the upper two-component spinor in fy G . In the chiral representation of the Dirac 
matrices, we require s 2 (r) +p 2 (r) > so that h\ 2 is invertible. As mentioned in section V 
this is a restriction on our variational ansatze. Using the known matrix elements for the 
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spin-isospin operators like T{X^ we then project eq. (jB6|) onto grand spin channels and obtain 
the desired second order differential equations. They may be written in the form, 

4 f f) ~\ 

^2\D G (r) + N G (r)— + M G (r)\ 9j (r) = , (B7) 



with 

G = diag ((G + 1)(G + 2), G(G + 1), G(G + 1), (G - 1)G) , (B8) 

and k 2 = uo 2 — m^. The matrices N G (r) and Mc(r) are given in terms of the functions 
s(r), . . . , 7(r) that specify the static background fields in the spherical ansatz, as in eq. (|2*7j) . 
Their elements are rather lengthy and we refrain from explicitly displaying them here. As 
r — > oo, N G (r) — > and Mc(r) — > and the differential equations decouple, as long as the 
potential goes to zero sufficiently fast. 

We have a four-channel scattering problem. We express the four wavefunctions and 
four boundary conditions in matrix form, Qij{r) = g^\r), where the linearly independent 
boundary conditions are labeled by j = 1,2,3,4. We then write Q{r) as a multiplicative 
modification of the matrix solution to the free differential equations, Q(r) = F(r) ■ H(kr), 
where H(x) = diag(/i^| 1 (a;), h^{x), h G (x), /i^i-^x)) with hf\x) denoting spherical Hankel 
functions of the first kind such that D G (r) -H(kr) = 0. (The 4x4 matrices F and H depend 
on the grand spin quantum number G. For convenience we omit that label from now on.) 
Imposing the boundary conditions F(r — > oo) = and F'(r — > oo) = 0, it is clear that the 
i th row of Q describes an outgoing spherical wave in the i th channel. Similarly, Q* describes 
incoming spherical waves. The scattering wavefunction can be written as 

G sc (r) = -g*(r) + g(r)S(k), (B9) 

and requiring this to be regular at the origin gives the scattering matrix 

S(k) = MmH-^k^F-^^F^^H^kr) . (BIO) 

r— >0 

We are interested in the sum of the eigenphase shifts in a given grand spin channel, 

5(k) = —TrlnS(k) = — limTr In (F^r)^*^)) • ( B11 ) 
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An efficient way to avoid any ambiguities in additive contributions of multiples of tt in S(k) 
is to define 

S(k,r) = — Trln (F _1 (r)F*(r)) , (B12) 



-3Tr (F'F- 1 ) (B13) 



2i 

with 5(h) = 5(k, 0). We then integrate 

d5(k, r 
Or 

along with F(k,r) from infinity to with the boundary condition ]im. r ^ 00 8(k,r) = to 
obtain S(k) as a smooth function of k. The differential equation for the matrix F(k, r), 

= F" + -F' + 2F'L' + \[F, O] + N(F' + FL') + MF , (B14) 

is obtained from 

d } 1 

D(r) + N(r)— + M(r) > Q(r) H~ 1 (kr) = Q, (B15) 
or J 

where L(kr) = In H(kr) and primes denote derivatives with respect to the radial coordinate. 
The components of L'(kr) can be expressed as simple rational functions, which avoids any 
numerical instability that would be caused by the oscillating Hankel functions. 

To construct the Born series for S(k), we introduce F^ n \k,r) where n labels the order in 
the background fields in an expansion around the zero-classical-energy configuration with 
/(°)(r) = in eq. (|32jl . We obtain the corresponding differential equations 

= jrW" + + 2F m> L > + J_r F (i) ) i + N (D L > + M (D 

= F^" + -F^' + 2F( 2 )'L' + O] + iV« (>«' + F^L'j 

+N®V + M«F« + , (B16) 

where the matrices and M'*' are obtained from iV and M by expanding to order i in the 
deviation of the background fields from the above described zero-classical-energy configura- 
tion. We integrate these differential equations with the boundary conditions F^'(k, oo) = 
and F^>'(k, oo) = and obtain 

^(Jfe) = -3tr (F (1) (M)) , 

5( 2 \k) = -3tr (f^(M)-^F (1) (M) 2 ) . (B17) 

We eliminate the quadratic divergence from the vacuum polarization energy by subtracting 
these from S(k) and adding them back in as renormalized first and second order Feynman 
diagrams. 
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There still remains the logarithmic divergence whose elimination would require third and 
fourth order Born subtractions. These become considerably more complicated, so instead we 
use the limiting function approach as described in The idea is to subtract only the local 
contributions to the third and fourth Born approximants to the phase shift by identifying 
them with the divergent contributions to the third and fourth order Feynman diagrams. 
To this end we formally manipulate these divergent Feynman diagrams. To extract the 
local contributions we set the external momenta to zero and then integrate over the energy 
and the two spatial angles of the loop momenta, k^, such that a (regularized) integral over 
k = \k\ is left. We write its integrand in the form as in eq. (|22j) . 



in the A = gauge and where = $ — v denotes the deviation of the Higgs field from its 
vev. 

Thus we numerically determine the bound state energies, the phase shifts, their Born 
series, and the limiting function. These are all ingredients in the expression for the fermion 
vacuum polarization energy in eq. ([22)1 . 




where 




(B18) 
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